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I  INTRODOCTI(»i 


The  objective  of  this  research  program  has  been  to  develop  the 
theoretical  models,  design  methodology,  and  technology  needed  for 
optimally  applying  near- field  electromagentic  sensor  arrays  to  non¬ 
destructive  evaluation  (NDE)  and  robot  control.  This  program  was  a 
collaborative  effort  by  SRI  International  and  Stanford  University, 
supported  by  separate  contracts.  SRI  has  studied  several  types  of 
electronically  scanned  arrays  composed  of  inductive  sensors ,  whereas 
Stanford'  has  placed  emphasis  on  the  use  of  capacitive  sensors. 

At  SRI,  we  have  found  that  small  printed-circuit  single-turn  loops 
exhibit  good  sensitivity  when  used  as  sensors. ^  This  technology  allows 
us  to  fabricate  high-resolution  arrays  readily.  By  addressing  different 
elements  in  the  array  and  suitably  processing  the  resulting  signals,  we 
can  realize  different  sensing  functions  with  the  same  array.  In  parti¬ 
cular,  we  have  demonstrated  the  capability  of  using  such  arrays  for  edge 
tracking  and  ranging  (proximity  sensing). 


l-l.tSJVtHJ 


^UTiT 


II  (CDELING 

An  inductive  sensor  array  located  over  an  imperfectly  conducting 
half-space  can  be  modeled  using  the  same  Fourier- transform  principles 
applied  by  Dodd  and  Deeds  to  model  eddy-current  coils. If  the  sensor 
elements  are  further  than  a  skin  depth  from  the  imperfectly  conducting 
surface,  it  is  also  possible  to  use  a  simpler,  but  approximate,  image 
model. We  have  developed  both  types  of  models  for  rectangular  coils 
under  the  assumption  that  the  fields  produced  by  a  system  of  coils  can  be 
computed  by  superimposing  the  fields  produced  by  each  coil  in  the  absence 
of  the  others. 

A.  Fourier -Transform  Method 

1 .  Fields  Above  an  Imperfectly  Conducting  Half  Space 

Consider  a  current  source  like  the  rectangular,  single-turn  coil 
shown  in  Figure  1,  but  of  arbitrary  shape  with  a  current  density  vector 
i(f).  The  quasi-static  magnetic  vector  potential  in  free  space  for  such 
a  source  located  in  a  volume  v'  is  Ao(f )exp( jwt)  where 


Aq  (r) 


i(r') 


and  /jg  is  the  magnetic  permeability  of  free  space. 

If  the  half  space  along  the  negative  z-axis  in  Figure  1  is  filled  with 
a  homogeneous  material  of  conductivity  a,  the  new  quasi-static  vector 
potential  A(f)exp(jwt)  can  be  found  by  solving  the  equation 

V^ACf)  -  -  +  jw/io(7A(f)  (2) 


FIGURE  1  RECTANGULAR  SINGLE-TURN  COIL  PARALLEL  TO  SURFACE  OF  CONDUCTING 
HALF  SPACE 

subject  Co  the  boundary  condition  that  the  x  and  y  components  of  the 
electric  and  magnetic  fields  be  continuous  at  z  -  0.  The  solution  to  Eq 
2  is  exponential  in  z.  The  x  and  y  dependencies  are  usually  found  by 
performing  a  two-dimensional  Fourier  transform  to  k  space,  and  then 
inverse  transforming  the  solution  to  real  space.  The  two-dimensional 
Fourier  transform  has  the  form; 


a. (k, z) 


1_ 

27r 


Aj(p,z)  d^p 


i  =•  X ,  y 


(3) 


where  p  and  k  are  the  vector  components  of  the  position  vector  r  and  the 
wave  vector  in  the  xy  plane,  respectively  (see  Figure  1). 


From  Eq.  2,  the  following  equations  for  aj  are  obtained.  For  z  >  0: 


-k2ai  + 


J  ij  (p.z) 


and  for  z  <  0: 


■k2aj  +  — — 


-  J<^Po*i 


Following  Beissner  and  Sablik^ ,  Che  solution  to  Eq.  4  can  be  written  as 


af(k,z)  -  ai(k,z)  +  (k)e‘'‘* 


and  for  Eq.  5: 


af(k,z)  -  -Yj  (k)e'^^ 


In  these  equations 


A  -  y  k2  +  jwpo<7  .  k  -  |k| 


ij  (k,z)  -  i^(p',z')e‘j^'^'**^U'^'l  dv’ 


where  a,-(k,z)  is  obtained  from  a  two-dimensional  Fourier  transform  of 
Aq (r)  given  in  Eq.  1. 

The  unknown  coefficients  (k)  and  7j (k)  are  found  by  matching  tangen 
tial  components  of  the  electric  and  magnetic  fields  at  the  conducting 
surface  z  -  0.  From  this  condition,  Beissner  and  Sablik  derived  the 
following  expression  for  the  magnetic  vector  potential  in  the  conductor 


4 


(10 


Aj  (p.z)  “  2^  J  ej^'P  d^k 

where 

<  -  ^  r  -  - 


(11 


and 


a^(k,z)  -  -  ^  ax(k,z) 


aj(k,z)  -  0 


(12 


(13 


However,  in  analyzing  an  electromagnetic  sensor  array,  the  electric  and 
magnetic  fields  above  the  conducting  surface  are  of  particular  interest 
These  fields  can  be  determined  by  matching  the  expression  for  af(k,z) 
(Eq.  6)  and  af(k,z)  (Eqs.  11-13)  at  z  -  0.  The  resulting  equations  are 


a^(k,z)  -  aO(k,0)|e'‘*  +  — 


k  •  A 


k+A 


(14 


ay (k,z) 


^  a^(k,z) 


(15 


and 


a. (k,z)  -  0 


(16 


The  electric  and  magnetic  fields  above  the  conducting  surface,  and 


are  then  determined  from 


5 


where 


E  (r)  -  -  jw  A  (p,z) 


H^(r)  -  —  V  X  A^(p,z) 
^0 


Aj(p,z)  -  —  aj(k,z)e^^*^  d^k 

2  IT 


Now  consider  the  special  case  shown  in  Figure  1,  where  we  have  a  single 
turn,  infinitely  thin,  rectangular  coil  carrying  a  current  I  in  a  plane 
parallel  to  the  surface  z  -  0.  To  find  the  electric  and  magnetic  field 
for  this  configuration,  we  must  find  the  current  density  i^  so  that 
a°(k,z)  can  be  calculated  from  Eq.  9.  These  current-density  components 
can  be  written  as 

i^(x,y,z)  -  I5(z-Zo)|^5j^y-(yo-a)  -  5  j^y- (yo+a)j  j  rectj^^^^j  ,  (20) 


i  (x,y ,z) 


I5(z-Zo)^5j^x-(xo-b)  -  S  x-(xo+b)jj  rect|^^j 


i^ (x,y,z)  -  0 


Inserting  these  expressions  into  Eq .  9  we  obtain 


a5(k,z)  -  j  r —  sin(k„a)  sin(kyb)  e‘j^’^0  e'^'l^'^O 


(k, z)  -  0 


Combining  Eqs .  23  and  14  and  inserting  the  resulting  a,-  (k,z)  into  the 
integral  of  Eq.  19  results  in  very  complex  integrals  because  of  the 
factor  (k-A)/(k+A)  that  occurs  in  the  expressions  for  af(k,z).  This 
factor  also  appears  in  the  image  theory  for  this  problem^ ,  where  it  is 
shown  that  it  can  be  approximated  by 


when  the  receiver  (sensor  array)  is  more  than  a  skin  depth  away  from  the 
conductor  surface.  The  factor  d  in  this  expression  is  related  to  the 
skin  depth 

5  -  y  2/(xjtMQa)  (27; 

by  the  relation 

d  -  (l-j)6  .  (28; 

Using  the  approximation  given  by  Eq.  26  we  obtain 


ax(k,z)  -  2a5(k,0)  exp(-kd/2)  sinh[k(z  + 


For  a  perfect  conductor  (cr  -*  <»)  ,  d  -  0  and  Eq .  29  simplifies  to 


a^(k,z)  -  2a2(k,0)  sinh(kz) 


7 


Comparing  Eqs .  29  and  30,  we  note  that,  if  the  factor  exp(-kd/2)  is 
included  in  the  expression  for  aJCk.O)  in  Eq.  23,  the  approximate  solu¬ 
tion  for  the  case  of  a  conducting  half-space  with  finite  conductivity  is 
the  same  as  that  for  a  perfectly  conducting  half -space,  except  that  z  is 
replaced  by  z  +  d/2 .  This  result  will  be  discussed  further  in  a  sub¬ 
sequent  section  on  image  theory. 

2.  Fields  Above  a  Perfectly  Conducting  Half  Space 


When  ax(k,z)  in  Eq.  30  for  a  perfect  conductor  is  used  in  Eq.  19,  the 
following  expressions  for  the  electric  and  magnetic  fields  are  obtained 
from  Eqs.  17  and  18  for  z  >  0: 


Ex(x,y,z)  -  -  —  sinh(kz)  a2(k,0)ej^'P  d^k 


I 

Ey(x, y,z)  “  ~  ^  sinh(kz)  a2(k,0)ej^’^  d^k 


Ej(x,y,z)  -  0 


(x,y,z) 


— —  cosh(k2)  a9(k,0)ej^'P  d^k 


Hy (x,y,z) 


(x,y,z) 


— —  k  cosh(kz)  a?(k,0)ej^’^  d^k 


The  integrations  appearing  in  these  equations  are  evaluated  using 
sine  and  cosine  tranforms  that  are  discussed  in  Appendix  A.  The  results 


rS'rWlTr.  ’W'JI  -U( 


i  ^Mn  I 
“  ^  4jr 


Fx (-z.p; , -p;)  +  Fx (-2, -p; . -p;)  + 


Fx(z,p;,p^)  +  Fx(z,-p;,pp 


Fx(-Z,p;,p')  +  Fx(-z,-p^ ,p')  + 


Fx  (z  .Pi  . -fiy  )  +  Fx  (z,  -PjI 


where 


Fx (±2,±  Pi ,±pi)  -  in 


a  ±  Pi  +  y  (a  ±  p;)2  +  (b  ±  p^)^  +  (zq  ±  z)‘ 
y  (b  ±  p')2  +  (zq  ±  z)^ 


and  p^  and  p^  are  the  x  and  y  components  of  the  vector  p-po .  Ex 
consists  of  eight  terms  containing  all  possible  combinations  of  the 
variables  ±z ,  ±pi ,  and  fp^ .  Ey  is  obtained  by  changing  the  sign  of  the 
right  side  of  Eq.  37  and  making  the  following  substitutions  in  Eqs .  37 
and  38; 

a  -►  b,  b  -►  a,  pi  -*  Py  ,  and  p'y  -*  p'n  (39^ 

Finally , 


Ej  -  0 


Similarly,  Hy  can  be  written  as 


I  r  sign  sign  ] 

Hy  -  ^  [  I  Gy(±z,±pi,-p;)  -  I  Gy(±z.±p;.p;) 


where 


Gy  (Iz.tpjl  ,1/)^  - 


(zq  ±  z)  (a  ±  />;) 


']/ 


(b  ±  p')^  +  (zq  ±  z)'^  y  (a  ±  p;)^  +  (b  ±  p')^  +  (Z(; 


z)‘ 


(42) 


sign 

The  symbol  ^  in  Eq.  41  denotes  the  sum  over  all  possible  sign 
combinations  of  the  variables  in  Gy  so  that  each  sum  appearing  in  Eq .  41 
consists  of  four  terms.  is  obtained  by  making  the  substitution  shown 

in  Eq.  39  in  Eqs .  41  and  42.  Finally, 


-  0 


(43) 


B .  Image  Method 

1 .  Model 

The  Fourier-transform  method  used  in  the  previous  section  is  quite 
powerful  and  provides  an  exact  solution  to  the  problem  of  calculating 
the  fields  generated  by  a  coil  over  an  imperfectly  conducting  half 
space.  However,  this  method  leads  to  integrals  that  must  be  evaluated 
over  infinite  domains  in  the  complex  plane,  which  may  be  difficult  to 
do.  Therefore,  for  such  half -space  problems,  it  is  usually  simpler  to 
use  image  theory.  Of  course,  image  theory  is  only  exact  for  perfect 
conductors,  but,  as  has  been  mentioned  in  Section  II -A- 1,  it  can  be 
extended  to  apply  to  imperfect  conductors  in  an  approximate,  but  simple, 
way . 

To  develop  the  image  theory  for  a  rectangular  coil  lying  in  a  plane 
parallel  to  a  conducting  half  space,  we  again  refer  to  Figure  1.  The 
current,  I,  flowing  in  the  coil  produces  a  vector  potential.  A"*",  in  the 
absence  of  the  half  space  given  by 


The  quantities  ,  ty  in  Eq.  (45)  are  the  x  and  y  components  of  r' . 
The  integrals  in  Eqs .  44  are  easily  evaluated  under  the  quasistatic 
approximation  kR  «  1.  The  results  are 


■  FUl  r.i  .tif  '»,■  v*  V  'v*  U-*  V.'^V.*  vr'^-’^V'Titvvr..'’ 


■.  ■■^.  '•■J  “■.■»■.' V.’'7'J'SrT^% 


+  Mol 
A^(r)  -  ^  in 


and 


+  _  Mol 

~  UT 


[R(r'  ;-a,l 


-][' 


b)  -  (r;-a)  R(r';a,-b)  -  (rJ+a) 


R(r';a.b)  -  (r^+a)  R(r':-a,-b)  -  (r^-a) 


(46a) 


R(r';-a,-b)  -  (r;-b)  ! 


R(r';a,b)  -  (r 


;+b) 


R(r':-a.b)  -  (r;+b)  R(r';a.-b)  -  (r;-b) 


(46b) 


where 


R(r';±a,±b)  -  J  (r;  ±  aY  +  (r;  ±  b)-:  +  (z  -  zq)' 


(47) 


For  a  perfectly  conducting  halfspace,  the  vector  potential,  A',  due  to 
the  image  current  can  be  obtained  from  Eq.  46  by  letting  I  -  -I  and 
zo  -  -Zg.  The  total  vector  potential  is  just  the  sum  of  the  two  parts, 
viz .  , 


A(r)  -  A+(r)  +  A-(f)  (48) 

Finally,  the  electromagnetic  fields  can  be  found  from  the  vector 
potential  by  using  the  relations 

E(f)  =  -jwA(r)  .  (49a) 

and 

H(f)  -  (1/hq)  V  X  A(r)  (49b) 
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It  can  be  shown  that  these  equations  are  identical  to  those  obtained  in 
Section  II-A-2  using  the  Fourier- transform  method. 

For  an  imperfectly  conducting  half  space,  it  has  been  shown^  that  the 
fields  above  the  half  space  can  be  determined  to  good  approximation  by 
placing  the  image  at  a  complex  depth  -(zg+d)  where  d  is  given  by  Eq.  28, 
Similarly,  for  an  imperfectly  conducting  layer  of  thickness  t  on  top  of 
an  imperfectly  conducting  half  space,  an  approximate  image  model  can  be 
obtained  by  using^ 


d  -  (1  -  j)SQ  (50) 

where 

1  -  e-'it/d  -j  )5 

Q -  (51) 

1  +  e-^t/(1-j)6 

Although  approximate,  these  formulas  allow  us  to  model  imperfect 
half -space  and  layered  half-space  problems  in  a  very  simple  way. 

2 .  Induced  Voltage 

To  compute  the  voltage  induced  in  different  sensor  coils  by  a 
separate  exciting  (driver)  coil,  we  use  Faraday's  Law 


'Sense 


-jw 


(VxA)  •  dS 


(52) 


where  S  is  the  surface 
use  Stoke 's  Theorem  to 
conductor  that  defines 


area  enclosed  by  the  coil.  Alternatively,  we  can 
write  this  equation  as  a  line  integral  along  the 
the  coil ,  viz . , 
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This  latter  equation  is  simpler  to  evaluate  since  it  uses  A  directly. 

Note  that  throughout  our  development,  we  assume  that  wire  in  the  coil  has 
very  small  cross  sectional  dimensions  so  that  the  current  in  the  coil  can 
be  represented  by  a  delta  function  in  this  cross  section  plane. 

3.  An  Example  of  Using  a  Sensor  Array  to  Determine  Proximit 


The  response  of  an  inductive  sensor  array  depends  on  its  proximity 
to  (lift-off  from)  a  nearby  conducting  half-space.  Thus,  it  is  of 
interest  to  examine  the  possibility  of  inverting  the  array  response  to 
provide  a  quantitative  measure  of  the  distance  of  the  array  from  the  half 
space . 

As  a  simple  example,  consider  the  array  of  concentric  rectangular 
coils  shown  in  Figure  2.  The  outermost  coil  is  an  exciting  coil,  and  the 
inner  coil  is  a  sensor.  The  coil  axes  are  assumed  to  be  perpendicular  to 
a  perfectly  conducting  half  space,  and  the  coils  are  located  at  a 
distance  Zq  above  the  half  space  (if  the  half  space  is  imperfectly 
conducting,  the  approximations  discussed  in  Section  II-B-1  can  be  used  in 
analyzing  this  problem) .  The  problem  is  to  determine  zq ,  given  the  array 
response . 

From  Eq.  53,  the  voltage  induced  in  the  inner  coil  by  the  outer  coil 
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-aa.Zo)  -  A,,  (x' ,aa,  Zq  )j  dx' 


+  A^(aa.y',Zo)  '  A^(-“a.y' .Zq )  dy' 


+  image  terms 


The  vector  potential  produced  by  the  exciting  coil  is  given  by  Eq .  46. 
Substituting  the  values  of  and  A^  of  Eq.  46  into  Eq.  54  results  in  a 
nonstandard  integral  of  the  form 


2a 


EXCITING  COIL 


SENSOR  COIL 


FIGURE  2  ARRAY  OF  SQUARE  CONCENTRIC  SINGLE-TURN  COILS 


This  integration  has  been  carried  out  symbolically  using  a  computer 
program  called  MAXIMA,  and  the  result  is 

-  X  .  (56) 

To  avoid  a  need  for  system  calibration,  we  assume  there  are  two  sense 
coils,  one  inside  the  other,  and  calculate  the  ratio  of  the  voltages 
induced  in  them  by  the  exciting  coil.  For  this  calculation,  we  assume 
that  both  the  source  and  detector  impedances  are  large  enough  so  that 
back  reaction  of  the  sense  coils  on  the  exciting  coil  is  negligible,  as 
is  interaction  between  the  sense  coils. 

The  normalized  lift-off  distance,  Zq  =  Zq/b  is  shown  in  Figure  3 
plotted  versus  the  voltage  ratio  computed  using  Eq.  54.  The  scale  of  the 
inner  sense  coil  was  taken  to  be  a  -  0.2,  while  the  size  of  the  outer 
sense  coil  was  allowed  to  vary  through  the  values  of  a  shown  in  the 
figure.  We  see  that,  the  larger  the  outer  sense  coil,  the  better  the 
resolution  for  determing  Zq  .  For  a  -  0.8,  the  lift-off,  Zg ,  can  be 
determined  most  accurately  for  values  of  Zg  between  about  0.03a  and  0.3a. 


C .  Reciprocity  Integral 

So  far,  our  modeling  has  been  concerned  with  a  relatively  simple 
geometry,  namely,  a  half  space.  It  is  also  important  to  model  situations 
where  the  sensor  array  is  scanned  parallel  to  the  surface  of  the  half 
space  in  order  to  detect  or  track  discontinuities  in  this  half  space. 

The  reciprocity  method  described  by  Auld^  is  useful  for  this  purpose.  We 
illustrate  this  method  for  the  case  where  there  is  a  step  in  the  surface 
of  a  perfectly  conducting  half  space. 


To  compute  the  result  of  scanning  a  perfectly  balanced  differential 
sensor  pair  across  a  small  step  in  a  conducting  plane  requires  that  the  E 
and  H  fields  produced  by  the  two  coil  configurations  shown  in  Figure  4  be 
computed.  Using  these  fields,  the  step  response  of  the  sensor  pair  can 
then  be  computed  by  evaluating  the  following  reciprocity  integral  for 
each  sensor-pair  position  (xQ.yQ.Zo): 


Va 


Zb 


e(1  )  X  h(2)  -  e(2)  X  )  Y 


a^  dxdy 


(57) 


In  this  equation,  the  superscripts  refer  to  the  two  cases  defined  in 
Figure  4  and  a^  is  a  unit  vector  in  the  z  direction.  The  integral  is 
evaluated  on  the  plane  z  -  0. 


CASE  1 


DRIVE  COIL  EXCITED; 
TEST  PIECE  PERTURBED 


CASE  2 


SENSOR  COILS  EXCITED: 
TEST  PIECE  UNPERTURBED 


FIGURE  4  CASES  USED  IN  RECIPROCITY  MODEL  FOR  COMPUTING  STEP  RESPONSE 


For  a  perfect  conductor,  Eq.  57  can  be  simplified  by  making  use  of  the 
vector  identity 


(e  .  h]  ■  i.  -  E  ■  (h  X  S,]  .  H  ■  (s.  X  e) 


Since  a^  x  E  is  proportional  to  the  component  of  E  tangential  to  the  x-y 
plane,  (E^^)  x  -  a^  -  0  everywhere  on  the  plane  z  -  0.  Similarly, 

(E^^ ^  X  a^  also  equals  zero  on  this  plane  for  y  <  0.  Hence,  Eq. 

57  reduces  to 
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This  formula  is  exact.  However,  there  are  no  analytic  formulas  for 
calculating  the  electric  fields  in  case  1,  and  approximations  are  often 
used.  The  approximation  we  have  used  to  compute  the  results  given  in 
Section  III -A  is  to  assume  that  there  is  no  step  in  the  conducting  half 
space,  but,  rather  that  the  half  space  is  displaced  downward  everywhere 
from  z  -  0  by  the  step  height.  This  should  be  a  good  approximation  if 
the  step  height  is  small  compared  to  the  radius  of  the  drive  coil  used  in 
case  1 . 

Similar  simplifications  to  Eq.  57  can  also  be  achieved  when  the 
stepped  half  space  is  imperfectly  conducting.  For  example,  if  we  assume 
that  the  tangential  fields  E^  and  H(  are  approximately  related  by  a 
scalar  surface  impedance®,  ,  viz., 

a^  X  E  =  ZjHt  ,  (60) 
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III  EVALUATION  OF  IMDUCTIVE  ARRAYS 
FOR  NEAR-FIELD  EDGE  TRACKING  AND  RANGING 


In  previous  work  on  this  project,  SRI  demonstrated  that  small 
single-turn  printed  loops  can  be  used  as  sensors  with  sufficient 
sensitivity  to  be  useful  in  NDE  and  robotics  and  that  printed-circuit 
techniques  facilitate  the  fabrication  of  arrays  of  small  loops  to  provide 
electronic  scanning  with  high  spatial  resolution.  To  further  evaluate 
this  capability,  we  have  built  such  an  array  and  have  evaluated  its  use 
for  closed-loop  edge  tracking  and  for  ranging  (proximity  sensing) . 

The  sensor  elements  used  in  this  work  consisted  of  single- turn 
rectangular  loops  with  inner  dimensions  of  0.030  by  0.050  in.  The  loops 
are  printed  on  circuit  boards  using  standard  techniques,  and  the  width  of 
each  printed  line  is  0.005  in.  A  typical  arrangement  of  two  rows,  each 
with  four  loops,  is  shown  in  Figure  5(a).  The  center- to-center  spacing 
of  adjacent  loops  in  a  row  is  0.070  in.,  and  the  distance  between  ad¬ 
jacent  printed  lines  of  two  neighboring  loops  is  0.010  in.  To  complete 
the  probe  system,  a  magnetic  drive  field  is  provided  by  a  30-turn,  0.5-in 
square  coil  (not  shown),  with  its  axis  perpendicular  to  the  circuit  board 
This  arrangement  results  in  a  two-port  probe  (called  a  reflection  probe 
in  eddy-current  testing);  the  drive  coil  is  the  probe's  input,  and  the 
sensor  loop(s)  is  (are)  its  output. 

To  detect  an  edge  or  groove  in  a  metal  surface,  it  is  best  to  mini¬ 
mize  lift-off  sensitivity  by  connecting  an  adjacent  pair  of  sensor  loops 
differentially  as  shown  in  Figure  5(b).  For  the  edge  -  tracking  experiment 
described  below,  we  used  the  signals  from  the  two  staggered  pairs  of 
sensor  loops  designated  as  Channels  1  and  2  in  Figure  5(a).  On  the  other 
hand,  for  ranging  (lift-off)  measurements,  we  used  the  signals  obtained 
separately  from  each  of  the  four  loops  in  a  row. 
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FIGURE  5  SENSOR  ARRAY 


Addressing  the  loops  in  pairs  or  individually  was  accomplished  using  a 
thin-film,  four-channel  read/write  amplifier  made  by  Silicon  Systems  (SSI 
520) .  Channel  selection  in  this  amplifier  is  controlled  by  the  data 
acquisition  computer.  The  two  1-kO  potentiometers  shown  in  Figure  5(b) 
are  used  to  balance  the  inputs  from  the  two  loops  (for  the  differential 
connection)  and  to  provide  the  proper  input  impedance  for  the  amplifier. 
For  single- loop  measurements,  the  resistance  of  each  potentiometer  was 


set  to  100  Q. 


A  block  diagram  of  the  measurement  system  used  on  our  experiments 
is  shown  in  Figure  6.  The  probe’s  drive  coil  is  excited  by  a  stable 
generator  operating  at  250  kHz  and  an  available  power  level  of  0  dBm. 

This  drive  level  gives  a  sensor-loop  voltage  level  of  the  order  of  10 
fiV .  The  output  from  the  differential  amplifier  is  filtered  by  a  250-Hz 
crystal  filter,  amplified  by  60  dB,  and  measured  by  a  Solartron  computing 
voltmeter  (Model  7151) .  The  level  of  the  signal  voltage  measured  by  the 
voltmeter  is  of  the  order  of  1  V,  and  the  corresponding  background  noise 


mSi 


level  is  about  1  mV.  The  digital  value  of  the  voltage  measured  by  the 
voltmeter  is  then  fed  to  an  HPIOOO  computer.  This  same  computer  also 
selects  the  channel  to  be  measured,  performs  any  required  computations  on 
the  measured  voltages,  and  using  a  suitable  algorithmic  interpretation  of 
these  voltages  as  a  basis,  controls  the  position  of  the  X-Y  stage. 

A.  Edge  Tracking 

Tracking  an  edge  requires  that  the  sensor  array  generate  a  discrimi¬ 
nation  function  such  that  a  change  in  the  sign  of  the  measured  voltage 
indicates  the  proximity  of  an  edge.  Such  a  discrimination  function  can 
be  generated  by  subtracting  the  outputs  of  Channels  1  and  2  shown  in 
Figure  5(a).  Three  discrimination  functions  measured  using  our  array 
are  shown  in  Figure  7 .  These  functions  were  obtained  by  scanning  the 
array  over  a  0 . 25- in. -deep ,  0.5-in. -wide  slot  in  a  flat  aluminum  plate 
(only  the  response  from  one  edge  of  the  slot  is  shown) . 

The  discriminator  voltage  is  plotted  in  Figure  7  as  a  function  of  the 
relative  position  of  the  array  center,  yo ,  with  respect  to  the  edge.  The 
array  center  was  located  as  close  as  possible  to  the  center  of  the  drive 
coil.  Because  the  two  sensor  pairs  that  make  up  Channels  1  and  2  are 
displaced  assyraetrically  with  respect  to  the  array  center,  the  discrimi¬ 
nation  function  crosses  zero  when  the  array  center  is  to  the  right  of  the 
edge  (about  0.035  in.  when  the  edge  is  perpendicular  to  the  scan  direc¬ 
tion)  .  The  slope  of  the  discrimination  function  depends  on  the  relative 
displacement  of  the  two  sensor  pairs,  which  is  0.070  in.  The  noticeable 
assymetry  of  these  functions  around  the  zero  crossing  is  caused  by  the 
nonuniformity  of  the  drive  field,  the  relative  setting  of  the  potentio¬ 
meters,  and  the  influence  of  the  other  edge  of  the  slot  at  large  values 
of  yo  . 


As  can  be  seen  in  Figure  7,  the  shape  of  the  discrimination  function 
varies  with  the  angle  0  between  the  scan  direction  and  the  edge.  When  0 
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OVER  METAL 


OVER  SLOT 


/  "  ^ 

i4 


SCAN  DIRECTION 


/ 


'  ^ 

i  V  / 

\  A  / 

\  \  \  / 

'  / 

1  \  /  / 

1  \l  /^60’ 


-200  -100  0  100  200  300  400 

RELATIVE  ARRAY  POSITION,  Vq —  thouiandth*  of  an  inch 

FIGURE  7  MEASURED  DISCRIMINATOR  CHARACTERISTICS  FOR 
DIFFERENT  EDGE  ROTATIONS  (CHANNEL  1  - 
CHANNEL  2) 


decreases,  the  two  peaks  of  this  function  move  apart  by  an  amount  propor¬ 
tional  to  1/sin  0,  and  the  peak  values  decrease.  This  kind  of  behavior 

1 

occurs  because  the  connections  that  are  used  in  these  experiments  produce 
an  array  that  is  polarized  along  the  rows  of  sensor  loops. 

A  simple  edge  -  tracking  algorithm  was  programmed  into  the  computer  to 
demonstrate  experimentally  that  such  a  sensor  array  can  be  made  to  follow 
an  edge.  The  resulting  measured  edge- tracking  characteristics  of  the 
system  are  shown  in  Figures  8(a)  and  8(b)  for  edge  angles  of  0  =  60“  and 
30°,  respectively.  These  characteristics  were  produced  as  follows;  the 
search  for  the  edge  starts  with  a  sensor  location  corresponding  to  the 
left  side  of  the  zero  crossing  in  the  discrimination  curve.  The  sensor 
steps  along  the  horizontal  y  axis  (0.010- in.  steps  were  used  in  this 
case)  until  the  discriminator  voltage  changes  sign.  At  this  point,  the 
sensor  steps  in  the  positive  x  (vertical)  direction.  If  this  step 
increases  the  discriminator  voltage,  the  scan  direction  is  reversed; 
otherwise  the  scan  continues  in  the  same  direction.  The  x  scan  continues 
until  the  discriminator  voltage  changes  sign  again.  At  this  point,  the 
probe  steps  one  step  in  the  y  direction  and  the  cycle  is  repeated. 

The  staircase-shaped  curves  shown  in  Figure  8  depict  the  endpoints  of 
the  probe's  movement;  any  backtracking  that  is  required  to  reach  these 
endpoints  is  not  explicitly  shown.  There  was  no  backtracking  in  the  step 
sequence  depicted  in  Figure  8(a).  For  this  case,  the  probe  stepped 
vertically  along  the  discriminator  curve  from  A-*A'->B,  then  horizontally 
from  B-*C ,  then  vertically  from  C-*D,  and  so  forth.  On  the  other  hand,  the 
step  sequence  shown  in  Figure  8(b)  did  involve  backtracking.  One  verti¬ 
cal  step  was  required  to  go  from  A-»B.  Then  one  horizontal  step  moved  the 
operating  point  along  the  discriminator  curve  from  B-*C  with  no  change  in 
sign.  However,  the  next  vertical  step  moves  one  from  C-»C' ' ,  which  is  in 
the  wrong  direction.  The  probe,  therefore,  reverses  direction,  moving 
the  operating  point  from  C'-»C' ' --because  of  backlash,  C' '  and  C  do  not 
coincide.  Because  no  change  in  sign  occurred  as  a  result  of  this  last 
step,  another  vertical  step  is  taken,  which  results  in  the  move  C''-*D. 
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Both  of  these  characteristics  were  generated  in  real  time,  thereby 
demonstrating  the  feasibility  of  tracking  an  edge  using  such  a  sensor 
array.  No  significant  difference  in  tracking  was  noted  between  the  30° 
and  60°  cases,  even  though  the  respective  discrimination  functions  are 
different.  Similar  tracking  characteristics  were  also  generated  using 
0.002-in.  steps,  thereby  demonstrating  that  excellent  signal - to-noise 
ratios  were  obtained.  Finally,  the  straight  lines  drawn  in  Figure  8  are 
for  angle-reference  purposes  only;  they  do  not  represent  the  position  of 
the  edge.  Determining  the  edge  position  with  respect  to  the  track 
requires  additional  calibration  of  the  probe. 

In  addition  to  this  experimental  work,  the  step  response  of  the  array 
was  modeled  using  the  approximate  theory  described  in  Sectio  II-C.  The 
computed  differential  response  of  a  pair  of  rectangular  sensor  loops  to  a 
0 . 004- in. -  step  in  an  alluminum  plate  is  compared  with  the  corresponding 
experimental  results  in  Figure  9.  The  array  position  has  been  normalized 
to  one  half  the  inner  dimension  of  the  drive  coil.  In  this  case,  the 
peak  response  occurs  when  the  sensor  pair  is  located  symmetrically  over 
the  edge  of  the  step.  Because  the  measurement  is  uncalibrated,  the 
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experimental  and  theoretical  curves  were  arbitrarily  normalized  to  make 
the  peak  values  equal  for  this  comparison.  We  see  that  the  theoretical 
response  curve  is  noticeably  wider  than  the  measured  curve.  Determining 
the  reason  for  this  discrepancy  will  require  further  study. 


Experiments  were  also  performed  to  demonstrate  that  the  distance  from 
the  sensor  array  to  the  workpiece  surface  can  be  determined  from  the  mea¬ 
sured  outputs  of  the  array  elements  without  moving  the  array.  This  is 
different  from  the  case  of  a  single  absolute  eddy-current  probe,  in  which 
the  distance  of  the  probe  from  the  underlying  surface  can  only  be  deter¬ 
mined  if  the  probe  output  is  measured  at  several  different  distances 
above  the  surface.  Our  ability  to  determine  range,  using  a  stationary 
array,  is  based  on  the  fact  that  different  elements  in  the  array  sample 
the  spatial -frequency  content  of  the  drive  field  differently. 

To  eliminate  the  effect  of  the  measurement-system  transfer  function 
in  determining  range,  we  measured  the  ratio  of  two  voltages  sensed  by 
different  elements  in  the  array  as  a  function  of  lift-off  distance.  For 
maximum  range  resolution,  the  elements  used  to  generate  this  ratio  should 
be  located  in  regions  where  the  drive  field  differs  as  much  as  possible 
with  respect  to  spatial  variation.  We  used  one  row  of  four  loops  and 
formed  the  desired  ratio  by  dividing  the  sum  of  the  voltages  sensed  by 
the  two  loops  furthest  from  the  center  of  the  drive  coil,  by  the  sum  of 
the  voltages  sensed  by  the  two  loops  closest  to  the  drive  coil’s  center. 
The  resulting  measured  ranging  characteristic  is  shown  in  Figure  10,  in 
which  the  lift-off  distance  (normalized  to  one  half  of  the  inner  dimen¬ 
sion  of  the  drive  coil)  is  plotted  versus  voltage  ratio.  It  is  seen  that 
a  10-percent  change  in  voltage  ratio  corresponds  to  about  a  50-percent 
change  in  lift-off  distance  for  a  normalized  lift  off  of  less  than  about 
0.3.  This  sensitivity  is  usable,  but  not  large.  We  expect  that  a 
different  array  geometry  might  improve  this  sensitivity. 

Using  the  image  theory  discussed  in  Section  II-B,  we  computed  the 
theoretical  ranging  characteristic  for  this  array  and  also  plotted  it  in 


IV  SUMMARY 


It  has  been  demonstrated  that  an  array  of  small  single -turn  printed 
sensor  loops  can  be  used  for  high-resolution  edge  tracking.  Also,  with 
proper  calibration,  absolute  range  can  be  determined  from  a  ratio  of 
sensor  voltages  without  moving  the  array.  Other  possible  uses  for  such 
an  array  are  in  material  characterization  and  two-dimensional  eddy- 
current  imaging.  Finally,  an  approximate  model  for  computing  the  array 
response  has  been  developed,  but  further  work  to  improve  its  accuracy  is 
required . 
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APPENDIX  A 

TRANSFX)RMATION  FROK  k  SPACE  TO  REAL  SPACE 

Transformation  of  the  electric  and  magnetic  fields  given  by  Eqs .  31-36 
from  k  space  to  real  space  is  obtained  by  transforming  and  Hy ,  and 
noting  that  Ey  and  can  be  obtained  by  making  a  simple  substitution  in 
the  expressions  for  E^  and  Hy  and  a  sign  change  in  Ej^  . 

Using  Eqs.  23  and  31,  E^  can  be  written  as, 
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r— j—  sinh(kz)  sin(k^a)  sin(ky b)ej  ^ ' ^0  ) e ^0  dk^  dky 
K  K)( 


If  we  write  this  equation  as 


I 


f  ( kx  ,  ky  )  e  j  ^  ^  ^  ■  ^0  )  dkx  dky 


and  use 


k  -  y  k^  +  k^ 


we  note  that  f(k^,ky)  is  even  in  k^  and  odd  in  ky  .  Hence,  the  only 
part  of  expressed  in  terms  of  sine  and  cosine  functions  that 

contributes  to  the  integral  is  jcos(k^p^ )sin(ky Py )  where  and  Py  are 
the  X  and  y  components  of  the  vector  p-pg .  Eq .  Al  can  therefore  be 
rewritten  as 


jJa*  A*  .1 


.  W/inI 


00 

f^f 

k  k,  [ 

1  • 


g-k (zq  -z  ) 


=  0+z)j  si 


in(kjj  a)  cos  (k^  p^  ) 


sin(kyb)  sin(kyPy)  dkj,  dky 


for  0  <  2  <  Zq - 


Using  the  identities 


sin(k_a)  cos(k„p’) 


sin(kyb)  sin(kyPy) 


j  sin(k^(a-Px)  +  sin(kj,  (a+pj| ) 
~  cos(ky(b-Py)  -  cos(ky(b+py) 


whe  expression  for  splits  up  into  8  terms,  each  of  which  can  be 
expressed  in  the  form 


WptQ  I 


Sx  “  -  j 


-30 

fl (±z ,±p; )  |cos  ky (b 


ky  (b  +  Py  )  j  j  -ikx  dky 


where 


1  sign  . 

fl(±z,±p;)  =  7— j—  X  -  z)  sin  k^(a  ±  p; ) 


slg 

The  symbol  2.  indicates  that  f^  is  a  sum  of  all  possible  sign 
combinations  of  z  and  p^  in  the  expression  on  the  right-hand  side  of  the 
equation.  Integrating  over  ky  we  see  that  the  integral  is  a  cosine 
transform  that  has  the  form  and  solution*’  shown  below  for  0  <  z  <  Zq ■ 


W5-I 


>>  ta’V'K^-JV'.V.^,^. VLV.’*."  VA'i 

v*.** 


CD 

(x^+a^)'^/2  exp  -  ^(x^+a^)'’/^  cos(xy)  dx  <=  Kq  a(/9^+y^)^/^  (A9) 


where 


x  =  ky,  a  =  kx,  /S-Zq  ±z,  y  =  b±p^ 


and  Kq (  )  is  the  modified  Bessel  function  of  the  second  kind  of  order 


The  remaining  integral  is  a  sine  tranform  and  has  the  form  and 
solution*^ 


x"'^  K„(ax)  sin(xy)  dx  = 


,  -  1.  .  ,]  r[,  -  .  1,) 


2  A  (;2  -  A 


[2 4-^- A  2  -A-/i 


(AlO) 


where 


x  =  k^  ,  y  =  a 


1,  /j  ■=  0 


(b  ±  p^)2  +  (Zg  ±  z)‘ 


r(  )  is  the  gamma  function,  and  j  (  )  is  a  hypergeometric  function. 
With  the  values  of  A  and  fi  given  above,  the  hypergeometric  function 
simplifies  to*^ 


(Ain 


y  (b  ±  +  (zq  ±  z)2 


Using  the  fact  that  r(h)  —  J  •k  ,  the  final  form  of  Ejj  becomes  that  given 
in  Eqs.  37  and  38. 

Eqs .  23  and  24  show  that,  in  the  integration  to  obtain  Ey ,  the  role  of 
and  ky  will  be  interchanged  from  that  used  in  the  integration  for  E^ . 
Hence,  Ey  can  be  obtained  from  Ey  by  performing  the  substitution  given 
in  Eq.  39,  and  changing  the  sign. 

In  deriving  the  expression  for  the  magnetic  field,  it  is  convenient  to 
start  with  Eq.  35  for  Hy .  Using  manipulations  similar  to  those  used  for 
Ey  ,  the  integration  can  be  expressed  as  the  following  cosine  transform*'^ 
in  ky  ; 


exp(-^(x^  +  a2)V2)  cos(xy)  dx  =  a{y^  +  /3^ )  ^ 


(A12) 


where  the  parameters  are  the  same  as  those  given  in  connection  with  Eq . 

A9 ,  and  (  )  is  the  modified  Bessel  function  of  the  second  kind  of  order 
one.  The  integration  over  ky  involves  the  sine  transform*^ 


CO 

Kj^(ax)  sin(xy)  dx  =  -  esc (-i/w)  (a^  +  y^  ) '  ^ 

4  2 

0 


-I  (y2  +  ) '  /2  +  y 


(y^  +  a^  )  ' 


(All) 


where  i/  -  I ,  a  -  [  (b  ±  Py)^  +  (z  ±  and  y  -  a  ±  py  .  Using  this 
result,  the  y  component  of  the  magnetic  field  can  be  written  as  shown  in 
Eqs.  41  and  42.  Hy  is  obtained  using  the  same  substitution  that  is  used 
to  convert  the  equation  for  E,  to  that  for  E,,  . 
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